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ABSTRACT 


The  effect  of  noise  variations  on  two  of  the  short  period  seismic  discriminants, 
spectral  ratio  and  log -likelihood  ratio  is  studied.  The  mean  square  errors  of  two 
different  estimators  of  these  two  discriminants  are  evaluated  and  some  theoretical 
and  experimental  results  for  two  seismic  signals  in  seismic  noise  are  presented. 
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I.  INTRODUCTION 


During  the  last  ten  years,  a  considerable  amount  of  effort  has  been  directed  to¬ 
ward  the  study  of  natural  and  manmade  seismic  sources  and  their  differentiating  fea¬ 
tures.  Based  on  theoretical  and  experimental  investigations,  a  number  of  discriminating 
features  of  earthquakes  and  underground  nuclear  explosions  have  been  discovered,  and 

their  identification  capabilities,  on  the  basis  of  data  from  large  magnitude  events,  have 

1,2 

been  demonstrated  .  When  the  event  magnitude  diminishes,  the  level  of  the  back¬ 
ground  noise  at  receiver  becomes  significant  and  causes  random  variations  in  the 
seismic  discriminants.  Thus  a  discriminating  feature  which  may  be  capable  of  dif¬ 
ferentiating  between  large  magnitude  seismic  events,  becomes  sensitive  to  noise  varia¬ 
tions  at  small  signal -to -noise  ratios  and  its  identification  capability  declines  drastically. 
Therefore,  the  identification  capability  of  each  seismic  discriminant  must  be  studied 
as  a  function  of  signal -to-noise  ratio.  This  study  requires  the  evaluation  of  the  statis¬ 
tics  of  each  seismic  discriminant. 

In  this  report,  we  will  consider  two  of  the  short  period  seismic  discriminants, 
spectral  ratio  and  log-likelihood  ratio,  treat  them  as  random  variables  and  study  their 
statistics.  By  statistics  of  a  random  variable,  in  this  study,  we  mean  the  mean  and 
mean  square  error  of  that  variable. 

11.  DEFINITIONS  OF  SPECTRAL  RATIO  AND  LOG -LIKELIHOOD  RATIO 

Spectral  ratio  seismic  discriminant  is  defined  as  the  ratio  of  energy  in 

quency  band  1.  45-1. 95  Hz  to  that  in  the  frequency  band  0.  35-0.  85  Hz  of  the 

energy  spectrum  i.  e.  ,  7 

SR  =  -£*— 

"L 


the  fre- 
P-wave 

(1) 
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where  Z,  and  ZH  are  the  energy  in  the  low  and  high  frequency  bands.  This  definition 

2 

differs  from  the  one  given  by  Lacoss  in  that  he  uses  the  voltage  spectrum  instead  of 

energy  spectrum.  The  primary  reason  for  the  use  of  energy  spectrum  is  that  it  is 

mathematically  more  convenient  to  study  the  statistics  of  the  energy'  spectrum  estimates 

than  those  of  the  voltage  spectrum  estimates.  In  the  absence  of  noise,  the  situation 

which  is  approximated  by  very  large  events,  ZL  and  ZH  depend  on  the  event  magnitude, 

source  mechanism,  transmission  path  and  etc.  Since  SR  is  a  ratio, we  can  consider 

normalized  ZL  and  ZH.  To  compute  normalized  ZL  and  ZH  for  a  large  magnitude 

event,  we  normalize  that  event  by  its  maximum  amplitude  on  the  unfiltered  trace  and 

then  compute  its  energy  in  the  low  and  high  frequency  bands.  Figs.  1  and  2  show 

scatter  diagrams  of  log10  ZL  and  logi0  ZH  for  all  of  the  shallow  earthquakes  and  pre- 

2 

sumed  explosions  of  ^  5.0  used  by  Lacoss  .  The  straight  lines  in  each  figure  are 
the  least  square  fits  for  the  earthquake  and  presumed  explosion  populations.  If  we 
define 


=  logio  ZL 

(2) 

S  lOglo  ZH 

(3) 

and  assume  that  at  a  given  magnitude  RL  and  RH  are  random  variables,  then  the  mean 
and  standard  deviation  of  each  of  these  variables,  for  a  given  population,  are  obtained 
by  the  height  of  the  line  corresponding  to  that  population  at  the  given  magnitude  and  the 
population  scatter  about  that  line,  respectively. 

The  log-likelihood  ratio  of  these  two  variable  is  defined  by 

s  ln  p(R,  ,  Rh:  given  presumed  explosions)  /- 

P(Rl»  Rh  :  given  earthquakes)  ' 
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Fig.  1.  Scatter  diagram  of  the  energy  in  the  low  frequency  band 
for  42  earthquakes  and  28  presumed  explosions.  Each  event  has 
been  normalized  to  have  a  unity  maximum  amplitude. 
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Fig.  2.  Scatter  diagram  of  the  energy  in  the  high  frequency  band 
for  42  earthquakes  and  28  presumed  explosions.  Each  event  has 
been  normalized  to  have  a  unity  maximum  amplitude. 
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where  p(RL ,  RH :  given  a  population)  denotes  the  joint  density  of  RL  and  RH  for  that 
population.  If  we  assume  that  RL  and  RH  are  independent  normally  distributed  random 
variables,  then 

p(RL,RH:  given  presumed  explosions)  = 


1 

2ttctx  i_  ct  x  h 


(5) 


where  mXL,  mXH,  aXH  and  a%i  denote  the  means  and  standard  deviations  of  R,  and  RH 
for  the  presumed  explosions'  population.  The  joint  probability  density  of  RL  and  RH, 
given  earthquakes,  is  obtained  by  replacing  X  by  Q  in  Eq.  (5).  Substituting  Eq.  (5) 
and  its  dual  for  the  earthquakes  into  Eq.  (4)  and  expanding  it,  we  obtain 

Y  =  (R,  ~mQI)2  .  (R,  -mx,)2  +  (Ru  -  m0H)3 

ql  2  c2  x  L  2^^  q  h 


(RH  -  mX  h)2 

2^2xh 


+  in  2s  LiZon 
CTxl  ^XH 


(6) 


We  now  observe  that  if  ctxl  -  a0L  =  aL  >  aXH  =  aXH  =  ctqh  =  aH  and  mXL  -  m0H 
maL  ~  l.  =  i  then  v  reduces  to 

oT2 

rh  "  rl  +  constant 

which  is  the  same  as  the  log  of  the  spectral  ratio  except  for  a  constant.  This  result 
implies  that  the  spectral  ratio  is  a  simplified  version  of  the  log-likelihood  ratio  dis- 

3 

criminant.  Kelly"  has  applied  the  log-likelihood  ratio  as  a  seismic  discriminant  and 
he  has  concluded  that  this  discriminant  pie rforms  the  same  as  or  better  than  the  spectral 
ratio. 
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In  order  to  classify  a  new  large  magnitude  event  using  either  of  these  two  discrimi¬ 
nants,  we  must  measure  m^,  ZL  and  ZH ,  and  evaluate  the  corresponding  discriminant 
and  compare  it  with  a  threshold.  Note  that  the  threshold  for  the  spectral  ratio  as  well 
as  mXL ,  mXH,  m5L  and  m^  for  the  log -likelihood  ratio  are  functions  of  m^. 

Thus  far  we  have  assumed  that  the  noise  contributions  to  ZL  and  ZH  are  small  by 
considering  large  magnitude  events.  When  the  relative  noise  level  becomes  significant 
we  must  replace  ZL  and  ZH  by  their  estimates  and  use  the  estimates  of  SR  and  y  for 
classification  purposes.  Estimates  of  SR  and  y,  for  a  seismic  event  observed  in  noise, 
are  random  variables  whose  statistics  vary  with  signal-to-noise  ratio.  Evaluation  of 
the  statistics  of  these  two  discriminants  is  the  topic  of  the  next  section. 

III.  MEAN  AND  MEAN  SQUARE  ERROR  OF  SR  AND  y  ESTIMATES 

A  A 

In  order  to  evaluate  the  mean  and  mean  square  error  of  SR  and  y ,  (estimate  of  SR 
and  y)  we  make  the  simplifying  assumption  that  the  Fourier  coefficients  of  the  observed 
time  series, 

Xj  =  st  +  nj  ;  i  =  1,2, - ,M  (7) 

are  statistically  independent  complex  random  variables.  This  assumption  is  not  un¬ 
reasonable  because,  according  to  the  Spectral  Decomposition  Theorem,  when  the  time 

window  becomes  large  the  Fourier  coefficients,  Xj's,  become  independent  of  one  an- 

/\  /\ 

other.  Under  this  assumption  ZL  and  ZH  (estimates  of  ZL  and  ZH),  which  consist  of  the 
sum  of  | S i  |2  (estimate  of  |Sj  |2)  in  the  low  and  high  frequency  bands,  become  indepen¬ 
dent.  Consequently,  the  joint  probability  density  of  ZL  and  ZH ,  which  is  required  for 

A  a 

the  evaluation  of  the  statistice  SR  and  y  reduces  to  the  product  of  their  individual 
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densities.  In  other  words,  the  problem  of  evaluating  the  mean  and  mean  square  error 

A  /\ 

of  SR  and  y  reduces  to  the  problem  of  evaluating  the  probability  density  of  an  energy 
estimate  in  a  frequency  band.  A  detailed  discussion  of  this  problem  and  evaluation  of 
the  mean  and  mean  square  error  of  any  estimate  of  SR  and  y  is  given  in  the  author's 
doctoral  dissertation*^. 

In  the  next  section,  we  will  present  some  theoretical  and  experimental  results 
for  the  following  two  energy  estimators  (defined  in  Reference  4): 

(i) 

Conventional  energy  estimator  =  E  |Xj  |2 

and  (ii) 

i  i  alx  I2  2 

ALPI1A-C  energy  estimator  =  E[ |Xt |  -  cctj  exp  (-— 1 — )  ] 

where  c  and  a  are  arbitrary  constants  and  2ctj2  is  the  noise  power  at  the  i'th  frequency 
component.  Note  that  ALPHA-C  estimator  is  a  function  of  the  observed  value  as  well 
as  the  noise  statistics.  For  the  evaluation  of  numerical  results  in  this  report  we  have 
chosen  c  =  1.  45  and  a  =  0. 1. 

IV.  NUMERICAL  RESULTS 

/s  A 

For  the  numerical  evaluation  of  the  mean  and  mean  square  error  of  SR  and  y  as 
function  of  SNR,  we  need  the  following  quantities: 

(i)  the  noise  power  density  spectrum  which  is  chosen  to  be  the  estimate  of  the  power 
density  spectrum  of  500  seconds  of  a  seismic  noise  process  observed  at  a  short  period 
seismometer, 

(ii)  the  signal  energy  density  spectrum  which  is  chosen  to  be  the  energy  spectrum  of 
10  seconds  of  a  large  magnitude  seismic  signal  recorded  by  a  short  period  instrument. 
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and 


(iii)  the  frequency  bands  of  interest,  the  frequency  increment,  and  hence  the  number 
of  frequency  components  in  each  frequency  band. 

The  short  period  seismic  data  is  sampled  every  0.  05  seconds  and  a  frequency 
increment  of  10/1024  Hz  is  chosen.  The  low  and  high  frequency  bands  of  interest  are 
0.  35-0.  85  Hz  and  1.  45-1.  95  Hz,  respectively.  From  these  the  number  of  Fourier 
coefficients  in  each  frequency  band  is  obtained. 

Fig.  3  shows  the  noise  power  density  spectrum  as  well  as  the  energy  density  spectra 
of  LASA  short  period  beams  of  a  large  magnitude  presumed  explosion  (signal  no.  1) 
and  a  large  magnitude  earthquake  (signal  no.  2).  These  spectra  are  normalized  to 
have  unity  energy  in  the  frequency  band  of  0.  0-2.  0  Hz.  The  parameters  associated 
with  these  two  signals  are  as  follows: 

(i;  for  signal  no.  1  (presumed  explosion) 


m,  =  5. 75 
b 

SR  =  0.  732 

y  =5.2 

mx  L  =3.  28 

mqi  =4.  01 

mx  H  =3. 13 

mQ  H 

=  2.68 

o2xl  =  0.027 

o3  Q  L  =  0.  1 

a2  x  h  =  0.  027 

^2 

G  Q  H 

=  0.  123 

and  (ii)  for  signal 

no.  2  (earthquake) 

m,  =  6.  0 
b 

SR  =  0.  0246 

y  =  -16.  5 

mx L  =  3.  37 

mq  L  =4.  03 

mx H  _  3. 1 

mQH 

=  2.  68 

as  xl  =  0.  027 

ct%l  =0-1 

a2x  H  = 0.027 

S2qh 

=  0. 123 

A  A  /N 

A 

We  now  present  the  mean  square  errors  of  ZL  and  ZH ,  SR  and  y  as  functions  of 
SNR.  The  experimental  results  are  obtained  by  adding  the  same  signal  into  50  noise 
segments,  each  of  10  seconds  duration.  In  order  to  control  the  signal -to -noise  ratio, 
we  keep  the  noise  level  constant  and  vary  the  signal  amplitude. 
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Fig.  3.  Power  density  spectrum  of  500  seconds  of  a  seismic  noise  process 
and  energy  density  spectra  of  two  seismic  signals  recorded  by  a  short 
period  vertical  instrument.  Signal  no.  1  corresponds  to  a  presumed  explo¬ 
sion  and  signal  no.  2  to  an  earthquake.  The  total  energy  in  each  of  these 
spectra  is  unity. 
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In  the  process  of  evaluating  the  mean  square  error  of  SR  and  y  at  each  SNR,  we 
obtain  the  probability  density  of  the  energy  estimates  in  the  two  frequency  bands. 

A 

Fig.  4  shows  the  cumulative  probability  distribution,  P(Z),  of  the  conventional  and 
ALPHA-C  estimates  of  energy  in  the  low  frequency  band  for  signal  no.  1.  For  a 

A  /\ 

given  value  of  the  variable  Z,  say  Za,  P(Za)  gives  the  probability  that  Z^Za.  The 
actual  signal  energy  in  this  frequency  band  is  1.  58  and  the  signal  energy  to  noise 
energy  (SNR)  in  this  band  is  -2.  1  db.  In  this  figure  are  also  shown  the  histograms  of 
the  conventional  and  ALPHA-C  energy  estimates  obtained  from  experimental  results. 

The  discrepancy  between  the  theoretical  and  experimental  results  is  due  to  two 
factors  (i)  dependence  between  the  Fourier  coefficients  in  the  frequency  band  of  interest 
and  (ii)  the  difference  between  the  sample  moments  and  statistical  moments  of  a  ran¬ 
dom  variable. 

Another  by-product  of  this  investigation  is  the  evaluation  of  the  mean  and  vari¬ 
ance  of  energy  estimates  in  a  frequency  band.  Figs.  5  and  6  show  the  normalized 
(by  the  actual  value  of  energy)  root  mean  square  error  of  the  conventional  and  ALPHA-C 
energy  estimates  in  the  low  and  high  frequency  bands  for  signals  no.  1  and  2.  These 
figures  present  the  statistical  (theoretical)  as  well  as  the  arithmetic  (experimental) 
root  mean  square  errors  as  functions  of  the  signal  multiplier  which  is  closely  re¬ 
lated  to  the  SNR.  Based  on  these  results  we  conclude  that  the  performance  of  the 
ALPHA-C  energy  estimators  is  superior  to  that  of  the  conventional  estimator  at  all 
signal -to -noise  ratios. 

Figs.  7  and  8  represent  the  statistical  and  arithmetic  root  mean  square  error 
of  the  conventional  and  ALPHA-C  spectral  ratio  estimators  for  signals  no.  1  and  2. 

Based  on  the  results  of  Figs.  5,  6,  7  and  8,  we  observe  that  although  the  performance 
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Fig.  4.  Probability  functions  and  histograms  of  the  conventional  and 
ALPHA -C  estimators  of  energy  in  the  low  frequency  band  of  signal 
no.  1  (presumed  explosion)  at  a  signal-to-noise  ratio  of  -2.1  db. 
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Fig.  5.  Normalized  root  mean  square  error  of  the  conventional 
and  ALPHA -C  estimators  in  the  high  and  low  frequency  bands  for 
signal  no.  1  (presumed  explosion). 
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SIGNAL  MULTIPLIER 


Fig.  6.  Normalized  root  mean  square  error  of  the  conventional 
and  ALPHA -C  energy  estimators  in  the  high  and  low  frequency 
bands  for  signal  no.  2  (earthquake). 
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Fig.  7.  Root  mean  square  error  of  the  conventional  and  ALPHA -C 
estimators  of  the  spectral  ratio  for  signal  no.  1  (presumed  explosion). 
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Fig.  8.  Root  mean  square  error  of  the  conventional  and  ALPHA -C 
estimators  of  the  spectral  ratio  for  signal  no.  2  (earthquake). 
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of  the  ALPHA-C  energy  estimator  is  superior  to  that  of  the  conventional  energy  esti¬ 
mator,  the  ALPHA-C  spectral  ratio  estimator  does  not  perform  as  well  as  the  con¬ 
ventional  spectral  ratio  estimator  in  the  case  of  signal  no.  2.  This  is  perhaps  due 
to  one  of  the  properties  of  the  ratio  of  two  variables  that  may  remain  constant  while 
both  of  the  variables  are  changing.  This  is  a  poor  quality  of  spectral  ratio  as  a  seis¬ 
mic  discriminant. 

Figs.  9  and  10  show  the  statistical  and  arithmetic  root  mean  square  error  of  the 
ALPHA-C  and  conventional  estimators  of  the  log-likelihood  ratio  for  signals  no.  1  and 
2.  These  figures  indicate  that  the  performance  of  the  ALPHA-C  estimator  of  the  log- 
likelihood  ratio  is  superior  to  that  of  the  conventional  estimator  of  this  discriminant 
for  both  signals. 

Comparison  of  the  mean  square  error  of  the  spectral  ratio  estimates  and  the 

mean  square  error  of  the  log-likelihood  ratio  estimates  for  signals  no.  1,  2  and  a 

few  other  signals  indicates  that  the  log -likelihood  ratio  discriminant  reflects  the 
/\ 

variations  of  ZL  and  ZH  in  a  more  pronounced  manner  than  the  spectral  ratio  dis¬ 
criminant.  Based  on  this  property  of  the  log -likelihood  ratio  and  the  fact  that  its 
identification  capability,  for  large  magnitude  events,  is  the  same  as  or  better  than 
the  spectral  ratio  its  application  is  highly  recommended. 
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Fig.  9.  Root  mean  square  error  of  the  conventional  and  ALPHA-C  estimators 
of  the  log-likelihood  ratio  for  signal  no.  1  (presumed  explosion). 
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Fig.  10.  Root  mean  square  error  of  the  conventional  and  ALPHA-C  estimators 
of  the  log-likelihood  ratio  for  signal  no.  2  (earthquake). 
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